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We present a theory of electronic properties of HgTe quantum dot and propose a strain sensor based on a 
strain-driven transition from a HgTe quantum dot with inverted bandstructure and robust topologically 
protected quantum edge states to a normal state without edge states in the energy gap. The presence or 
absence of edge states leads to large on/off ratio of conductivity across the quantum dot, tunable by adjusting 
the number of conduction channels in the source- drain voltage window. The electronicj)roperties of a HgTe 
quantum dot as a function of size and applied strain are described using eight-band hp Luttinger and 
Bir-Pikus Hamiltonians, with surface states identified with chirality of Luttinger spinors and obtained 
through extensive numerical diagonalization of the Hamiltonian. 

There is currently significant interest in using quantum effects to develop capabilities in sensing at the 
nanoscale, from single-electron charge detection using quantum point contacts^'^, few nuclear spins using 
the quantum states of NV center in diamond^'^ to piezotronic sensors of strain^'^ and development of smart 
skin^. The piezotronic sensors rely on inducing a charge on the surface of piezoelectric semiconductor nanowire 
by strain, which in turn continuously changes the conductivity of the nanowire^. Here we describe a novel strain 
sensor based on a strain-driven transition from a HgTe quantum dot with inverted bandstructure and robust 
topologically protected quantum edge states to a normal state without edge states in the energy gap. The presence/ 
absence of edge states is expected to lead to large on/off ratio of conductivity across the quantum dot in analogy to 
recently demonstrated large piezoresistive current on/off ratios driven by a pressure-induced metal/insulator 
transition in rare earth chalcogenide thin fllms^ The electronic properties of a HgTe quantum dot as a function of 
applied strain are described using eight-band hp theory and surface states are identified with chirality of Luttinger 
spinors. 

The principle of operation of the quantum strain sensor based on a HgTe topological insulator quantum dot is 
summarized in Fig. L Panel (a) shows a thick, finite-size HgTe quantum well in the form of a disk embedded in 
vacuum or a higher bandgap material such as CdTe. The quantum dot is connected to two metallic electrodes. The 
quantum well of HgTe is an example of a topological insulator (TI), a material with an energy gap in the bulk, 
accompanied by helical, topologically protected states at its edge^"^^. This is in contrast with normal semi- 
conductor (NS) quantum dot, in which the states from the electron and hole subbands are separated by an empty 
energy gap^^. The thickness of the HgTe quantum well, through quantum confinenement, controls the relative 
position of the s and/> bands. As a result, there exists a transition from the inverted to normal bandstructure when 
the thickness of the quantum well is reduced below the critical thickness He. At this special thickness, the 
quasiparticle in-plane dispersion correponds to a single Dirac cone^l The thickness of the disk in Fig. 1 is larger 
than He, which results in the inverted bandstructure and the existence of states, whose energy falls within the gap 
and probability density peaks at the disk edge^^. 

In Fig. 1(a), the red ring represents the computed probability density of one of the edge states. If the Fermi 
energy of the leads is aligned with the energy of the edge state, electrons are expected to tunnel efficiently from the 
right to the left lead via the edge state and one detects a high current flow. Assuming ballistic transport, the current 
I is proportional to the number M of edge states and the applied voltage 1 = GqMTVsd^ where Gq is a quantum 
of conductance Gq = e^/h, and Tis the transmission coefiicient^^"^^. This is schematically visualized as the red line - 
high current or large number of edge states - in Fig. 1, where we plot the tunneling current as a function of the 
strain. As we apply strain to our disk, the edges of both s andp bands shift in energy, but with different deformation 
potentials. If the strain is large enough, the HgTe quantum well bandstructure is no longer inverted and the edge 
states disappear from the energy gap and are out of the resonance with the Fermi levels of the leads. Since there is 
no state available for the electrons to tunnel through, the current is expected to vanish, as depicted in Fig. 1 (b). As 
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Figure 1 | Schematic explanation of the principle of operation of the 
strain sensor, (a) The unstrained HgTe disk in the inverted regime is 
placed between two metaUic electrodes. The Fermi energy of the electrodes 
matches the energy of one of the edge states. The current flows since the 
electrons tunnel easily into and out of the edge state in the disk, (b) Upon 
applied vertical stress, the gap opens in the single-particle spectrum of the 
disk. There are no states matched with the Fermi energy of the leads, the 
current does not flow. 
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Figure 2 | (a) Subband edge energies in a HgTe quantum well Sit kx = ky= 0 
as a function of the well thickness H. He denotes the crossing between the 
lowest conduction (blue) and heavy-hole (red) subband edges, marking 
the transition between the normal semiconductor and the inverted 
bandgap material. Inset shows a schematic picture of the system composed 
of a single, free-standing, disk-shaped HgTe nanocrystal. (b) The in-plane 
dispersion of the conduction and heavy-hole subbands at the height 
showing the characteristic Dirac cone shape. 



shown in the main panel, the vanishing of current (blue line) is 
expected to occur when the strain exceeds the critical value. The 
robustness and sharpness of edge states is expected to lead to a very 
high on- off ratio of current as a function of applied strain. 

We now turn to a microscopic description of a HgTe quantum dot 
based on eight-band hp theory developed for semiconductor 
quantum dots^^'^°"^^ HgTe quantum wells^^ and HgTe colloidal 
nanocrystals^^"^^. We note that a simplified model of a strictly two- 
dimensional quantum disk embedded in vacuum and described by a 
heavy hole and conduction state has been recently studied by Chang 
and Lou^^. Specifically, we focus on a model colloidal HgTe quantum 
disk with radius R and height H in vacuum, shown schematically in 
the inset to Fig. 2(a). The single-particle energy levels and the cor- 
responding wave functions as a function of strain applied along the 

disk height are computed in the eight-band hp approach with the 
HgTe material parameters, described in detail in the Methods section 
and taken from Ref. 24. 

The eight-band hp bulk Hamiltonian is written in the basis of two 
conduction and six valence subband states. The eigenstates of the 
Hamiltonian are spinors with chirality up or down (see Methods 
section). Each spinor carries conduction band, heavy hole, light hole 
and spin split- off band components with a specific angular momenta 
for in-plane motion and parity for vertical motion. The chirality is a 
good quantum number and allows for the rigorous classification of 
quantum states as observed experimentally by Doty et aP^'^^. For 
example, we write the spinor with chirality "up" in the form 
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where, | n, m) denote the in-plane basis functions with radial node 
number n and angular momentum m, |2/ + 1) ( |2/)) describe the even 
(odd) trigonometric functions in the z-direction, |/, are the sub- 
band microscopic (Bloch) functions, and A^^^ are expansion coeffi- 
cients (see the Methods section). We see that the spinor contains all 
conduction band and valence band states and that the different sub- 
band components enter with different angular momentum and par- 
ity. This spinor is characterized by the total angular momentum 
3 

quantum number L = m+ - which is a sum of the orbital angular 

momentum m and the z-projection of the Bloch angular momentum 
/. We have chosen to define the angular momenta of all spinor 
components relative to the number m of the heavy-hole component. 
In a similar way we define chirality "down" spinors. The eight-band 
hp Hamiltonian is expanded in the basis of spinors and diagonalized 
to obtain coefficients A as explained in the Methods section. The 
wave functions vanish at the edge of the disk and any state localized 



SCIENTIFIC REPORTS | 4 : 4903 | DOI: 1 0.1 038/srep04903 



2 



at the edge must be characterized by its effective radius i^* and decay 
into the disk. Obtaining such locaUzed states in terms of our basis 
functions is a nontrivial numerical task. 

First we search for an optimal height H of the HgTe quantum well 
close to the transition from the normal semiconductor to the topo- 
logical insulator, i.e., from the normal to the inverted bandgap 
regime. Figure 2(a) shows the energies of the quantum well subbands 
in the vicinity of the critical thickness He, where the conduction band 
edge belonging to the lowest subband (denoted in blue) crosses the 
lowest- subband heavy-hole band edge (denoted in red). Thus, at 
sufficiently small well thicknesses H we deal with the normal semi- 
conductor phase with the positive bandgap, while at larger well thick- 
nesses the material exhibits a bulk-like, inverted bandstructure. The 
critical thickness for the infinite HgTe quantum well surrounded by 
vacuum is found at He = 3.4 nm. Here, the in-plane dispersion 
corresponds to a single Dirac cone as discussed and observed for 
HgTe/CdTe quantum welP^. This Dirac cone is presented in Fig. 2(b). 

We now turn to a quantum disk with a height of H = 4 nm which 
corresponds to a quantum well in the TI regime. We unfold our 
spinors in linear combinations of N Bessel functions and M sines 
for each subband, which results in the Hamiltonian matrix of order of 
K = S NM for each total angular momentum and chirality channel. 
In Fig. 3(a) we show the eigenenergies of the quantum disk obtained 
without strain by diagonalizing the Hamiltonian matrix with N = 40 
and M = 10. The levels, denoted by red and blue bars for chirality 
"up" and "down", respectively, are plotted against the total angular 
momentum L. The states form degenerate Kramers pairs, one with 
chirality "up", and the other with chirality "down", characterized by 
opposite total angular momenta. We see the formation of two bands 
of edge states with quasi-linear dispersion, the band with chirality 
"up" decreasing in energy as the total angular momentum is 
increased, and the band with chirality "down" behaving oppositely. 
Further, at energies higher than the edge bands we find the "interior" 
(non-edge) electron states, resembling those found in the NS 
quantum dot. The "interior" heavy-hole states form a ladder of levels 
below the edge bands. 

Fig. 3(b) illustrates the effect of strain included via the Bir-Pikus 
Hamiltoninan on the single-particle energies. The compression s = 
—0.02 along the disk height shifts differently the electron and hole 
levels and tunes the HgTe bandgap, such that its increase is analogous 
to moving from right to left in Fig. 2(a). This results in a wider energy 
gap in the quantum dot spectrum and removal of the edge states. 

We now discuss how the strain induced transition from a state 
with edge states in the energy gap to the insulating state without such 
states might be detected in transport. The two horizontal lines visible 
in Fig. 3(a) and (b) denote the conduction window AEsd = eVsn^ 
where Vsd is the source-drain voltage, defined by the Fermi levels of 
the left (L) and right (R)leads. The number M of edge states with the 
energy found within this window contributes to the tunneling of 
electrons from one lead to the other, and current I = GqMTVsd as 
discussed in the introduction. With our choice of the AEsd^ at zero 
strain we find six states, while at £ = — 0.02 the conduction window is 
empty and the tunneling current cannot flow. In Fig. 3(c) we show 
the effect of the width of the source-drain voltage window on the 
number M of available conduction channels as a function of the 
strain. Different curves correspond to different source -drain vol- 
tages. We see that a small conduction window contains fewer edge 
states, exhibits a low critical strain and an abrupt transition from the 
"on" to "off state. This is the quantitative illustration of the central 
result of our work presented in Fig. 1. As the conduction window 
increases, the "on" current increases due to the increased number of 
conduction channels, and the transition to the "off regime is less 
abrupt, consists of steps, and exhibits a larger critical strain. We 
expect, therefore, that the sensitivity and signal-to-noise ratio can 
be significantly tuned simply by changing the source -drain voltage of 
the leads. The effect discussed here is to be contrasted with piezo- 
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Figure 3 | Single-particle states of the quantum disk of height H = 4 nm 
and vertical strain e^^ = 0 (a), and Szz — ~0.02 (b). The black horizontal 
lines define the conduction window established by the Fermi levels of the 
leads of the electrical strain detector. Red (blue) bars correspond to 
chirality "up" ("down"). Panel (c) shows the current flowing through the 
sensor, measured in terms of the number of conduction channels, as a 
function of the strain. Different curves correspond to different widths of 
the conduction window. 

electric sensors where induced charge is a linear function of applied 
pressure^"^, but has analogies to single- electron charge detection 
using conductance steps in quantum point contacts ^'^ and pressure 
induced metal-insulator transitions in thin films^. The magnitude of 
strain of few percent is comparable to strains detected using piezo- 
electric semiconductor nanowires^. 

We conclude our proposal by demonstrating that not only the 
number, but also the nature of the conduction channels changes with 
the strain. In Fig. 4 we show the radial dependence of the probability 
densities of the state denoted in Fig. 3 as |A). The panels show, 
respectively, the electron (red) and heavy-hole (blue) component of 
this state. We see that as the strain is increased, the probability density 
of the electronic component evolves from one peaked at the edge to 
one peaked in the interior of the disk, while the hole component drops 
to zero. Thus, as expected, the states begin to resemble those of the NS 
quantum disk. As the probability density moves to the center, the 
tunneling from the leads onto the conduction channels decreases 
exponentially. Therefore, even though a state may be still within the 
conduction window, its contact to leads (tunneling matrix element), 
and therefore the tunneling current, decreases as a function of stress. 

In conclusion, we present a microscopic theory of a HgTe 
quantum dot as a function of size and applied strain. We demon- 
strated the existence of edge states and their removal with applied 
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strain. We propose that this mechanism could be used to implement 
a nanoscale, all -electrical, low-energy strain sensing device, in which 
the presence of a strain beyond a threshold value could be detected 
electrically as a collapse of the tunneling current through the edge 
states of the HgTe quantum disk. 

Methods 

The essential physics of the TI edge states can only be captured if one accounts for the 
strong mixing of the electron and hole subbands"'"^^'^*^"^". Here we analyze the single- 
particle properties of a HgTe disk, both free-standing and embedded in a normal 
semiconductor material in the eight-band hp approach. This allows to relate its 
electronic properties directly to its geometric and material properties. Building on our 
previous work with NS quantum wells^^ and dots^°"^^, we consider a HgTe quantum 
disk, shown schematically in the inset to Fig. 2(a). Here we discuss a free-standing 
nanocrystal, i.e., one whose surface is modeled simply as an infinite potential barrier. 
However, similar results were obtained for a HgTe quantum disk embedded in the 
CdTe barrier material. We take the radius of = 55 nm and vary the height H from 
2 nm to about 10 nm. As already mentioned, in our calculations we employ the eight- 
band hp approach. The Hamiltonian and relevant HgTe material parameters, 
described in detail in the following, are taken from Ref 24. We note that a strictly two- 
dimensional model of the HgTe dot has been considered in Ref. 16. In that approach, 
the effects of the disk height and subband mixing were accounted for only through 
effective parameters of the electron and the heavy hole in a two-band approach, which 
made its results difficult to relate directly to realistic structure parameters. 

The eight-band hp bulk Hamiltonian is written in the basis of two conduction and 
six valence subbands. Denoting the spin of the quasiparticle with an arrow, a = | (|) 
= ±1/2, the Bloch basis set is chosen in the following form: 

(2) 
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The total Bloch angular momentum for the electron subbands (the first two states) is 
equal to the spin and is 1/2. The total Bloch angular momentum for the hole states is 
3/2 for the heavy and light hole subbands, and 1/2 for the spin-orbit split-off sub- 
bands, and its projections are indicated by the second quantum number in the ket. 
The bulk Hamiltonian written in this basis takes the following form: 
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(10) 



S+ = —2\/3 y^k+'k^, 
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and k+ = kx± iky. Further, mo is the mass of a free electron and H is the Dirac 
constant. In what follows we take the following HgTe material parameters'"^: yi = 4. 1, 
72 = 0.5, 73 = 1.3, the spin-orbit splitting A = 1.08 eV, and the bandgap Eg = Ec — 
= —0.303 eV, appropriate for our negative-badgap material. The conduction-val- 
ence subband coupling parameter P can be deduced from the Kane energy 
Ep = 2mQP^ /H^ = 18.8eV. Lastly, we adopt the axial- symmetric approximation in 
which y = (72 + 73 ) /2 and the operator R can be written in a simpler form 

R = \/2,y k_ . The Hamiltonian ( 10) is applied to the zero -dimensional quantum 

2mo 

disk by the usual substitution kx= — id / dx, and analogously for coordinates y and z, 

and k+ = —ie-^'^ ± i-^-\ . The natural choice of basis for the single-particle 

\dr r d(p J 

states consists of eigenvectors of the single-band Hamiltonian, 
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where is the Bessel function of m-th order, and is the n-th zero of that function. 
The quantum numbers ?i = 0, 1, m = 0, ± 1, ±2, and / = 1, 2, ... are, 
respectively, the nodal number and angular momentum of the in-plane function and 
the vertical subband index. Note that even (odd) values of / correspond to odd (even) 
vertical wave functions, with the origin of the coordinate system placed in the center 
of the base of the disk. The functions by construction vanish on all the surfaces of the 
disk, and therefore working with this basis we do not have to enforce any additional 
boundary conditions. 

The eigenvectors of the Hamiltonian (10) are sought in the form of eight-com- 
ponent spinors, each component being a linear combination of single- subband 
functions. However, the analysis of symmetries of the Hamiltonian allows to discern 
two classes ("chiralities"), whose existence is related to the Kramers degeneracy. 

Owing to the fact that the operator k+(k-^ increases (decreases) the envelope 

angular momentum by one unit, and the operator kz flips the symmetry of the vertical 
subband, we write the spinor with chirality "up" in the form presented in Eq. 1 (Refs. 
20-22). This spinor is characterized by the total angular momentum quantum 
3 

number L = m H — which is a sum of the orbital angular momentum m and the z- 

projection of the Bloch angular momentum. We have chosen to define the angular 
momenta of all spinor components relative to the number m of the heavy-hole 
component. The spinor with chirality "down" takes the form 
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The operators appearing in the above matrix are defined as: 



where N and M define respectively the total number of radial and vertical harmonics 

3 

used in the calculation, and the total angular momentum L = m . The problem is 
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Figure 4 | Probability densities of the state | A) from Fig. 3 as a function of the radial coordinate for different applied strains. The electron and heavy 
hole wave function components are shown in red and blue, respectively, in the form of three-dimensional plots in the left-hand column. The radial plots 
along the x axis are shown on the right with the same colour code. Densities for the strain e = 0, 0.01, and 0.02 are shown in the top, middle, 
and bottom panels, respectively. 



thus reduced to finding the sets of coefficients A^^^^ and B^^^^ together with corres- 
ponding energy eigenvalues by diagonalizing numerically the Hamiltonian matrix set 
up in the above basis for each chirality, respectively. 

Following Novik et al.^'* we account for the strain effects by means of the Bir-Pikus 
Hamiltonian which is added to the eight-band Hamiltonian, Eq. (10). In the above 
basis, this Hamiltonian takes the following form: 
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with the following definitions of operators: t = ac (exx + %y + Szz) > 



U = ay {Bxx + %y + fizz) > V = - ^ (Sxx + % - 2£zz) ,S=-d{£xz- i^yz), 

a/3 

f = —h(&xx ~ %) + id&xy The strain defined by strain tensor matrix elements &ij 

is translated into energy via deformation potentials Uy, b, and d. In this work we 
utilize the values given by van de Walle^^: Uc = —4.60 eV, = —0.13 eV, and b = 
— 1.15 eV. The value for the potential d is not needed as we do not consider any 
nondiagonal (shear) strain in the system. 

In this work we consider a specific case of the disk being stressed in the vertical 
direction (along its thickness). We account for this case in the Bir-Pikus Hamiltonian, 
Eq. ( 18) by taking a nonzero strain tensor element b^z = AH/H, describing the relative 
change of the dot height, and setting all other strain tensor matrix elements to zero. 
This approach neglects any buckling effects that might occur close to the edges of the 
disk. Under such simple strain, the only remaining nonzero elements of the Bir-Pikus 
Hamiltonian are i, which renormalizes the conduction band edge, and u and v, which 
influence the valence subbands. Specifically, if the disk is stressed compressively (Szz 
< 0), the conduction band edge is shifted towards higher energies, since the relevant 
deformation potential < 0. The heavy-hole band edge undergoes the shift of 
w + V = (av — b)Exx. Both valence-band deformation potentials are negative, but the 
potential b is an order of magnitude larger than Uy. Under the negative strain, 
therefore, the heavy-hole subband edge is shifted towards lower energies. In total, the 
corrections to the conduction and heavy-hole band edges contribute a positive term 



SCIENTIFIC REPORTS | 4 : 4903 | DOI: 1 0.1 038/srep04903 



5 



to the bandgap, i.e., tend to drive the system from the inverted regime towards the 
normally gapped regime. This property is the physical principle of operation of our 
strain sensor. 
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